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Antiferromagnets (AFMs) exhibit intrinsic magnetization when the order parameter spatially 
varies. This intrinsic spin is present even at equilibrium and can be interpreted as a twisting of 
the homogeneous AFM into a state with a hnite spin. Because magnetic moments couple directly 
to external magnetic he Ids, the intrinsic magnetization can alter the dynamics of antiferromagnetic 
textures under such inhuence. Starting from the discrete Heisenberg model, we derive the continuum 
limit of the free energy of AFMs in the exchange approximation and explicitly rederive that the 
spatial variation of the antiferromagnetic order parameter is associated with an intrinsic magnetiza¬ 
tion density. We calculate the magnetization prohle of a domain wall and discuss how the intrinsic 
magnetization reacts to external forces. We show conclusively, both analytically and numerically, 
that a spatially inhomogeneous magnetic held can move and control the position of domain walls 
in AFMs. By comparing our model to a commonly used alternative parametrization procedure for 
the continuum helds, we show that the physical interpretations of these helds depend critically on 
the choice of parametrization procedure for the discrete-to-continuous transition. This can explain 
why a signihcant amount of recent studies of the dynamics of AFMs, including ehective models 
that describe the motion of antiferromagnetic domain walls, have neglected the intrinsic spin of the 
textured order parameter. 

PACS numbers: 75.50.Ee, 75.60.-d, 75.78.Fg 


I. INTRODUCTION 

Measuring the ordered state of antiferromagnets (AFMs) 
is complicated by the absence of macroscopic magnetiza¬ 
tion. The promise of AFMs as candidates for active roles 
in spintronics logic elements have increased the interest in 
addressing this problem.In particular, the observation of 
tunneling anisotropic magnetoresistance in AFMs^“^ repre¬ 
sents a clear experimental procedure to detect the antifer¬ 
romagnetic order. Furthermore, current-induced torques 
on the antiferromagnetic order have been theoretically pre- 
dicted^^^^ and experimentally indicated in spin valve sys¬ 
tems. Also the ferromagnetic concept of spin pumping 
has been generalized to AFMs.^^ The possibility of ma¬ 
nipulating the antiferromagnetic order parameter by exter¬ 
nal forces has fueled renewed theoretical interest in domain 
wall motion in AFMs due to both charge^^^^^ and spin^^“^^ 
currents. However, the reports on current-induced domain 
wall motion^^ are based on indirect observations and not 
confirmed by other methods or groups. Therefore, there is 
no straightforward method to reliably detect the dynamics 
of textures in the antiferromagnetic order. 

In this article, we discuss the intrinsic magnetization as¬ 
sociated with an inhomogeneous antiferromagnetic order 
parameter. We describe the origin of the intrinsic spin 
and discuss whether it can be exploited to detect antiferro¬ 
magnetic texture dynamics, e.g., domain wall motion. To 
revisit this topic, which was pioneered for one-dimensional 
systems in Refs. 20-22, we construct the continuum free 
energy functional for AFMs from the discrete Heisenberg 
Hamiltonian in the exchange approximation. We use the 
Hamiltonian approach to show that the intrinsic magneti¬ 
zation due to textures in the order parameter arises from a 
parity-breaking term in the energy functional that is ab¬ 
sent in a commonly used alternative parametrization of 
the continuum fields. We clarify the mapping between the 
two different parametrizations and explain how the intrin¬ 


sic magnetization can be easily missed in models which are 
based on the alternative parametrization. We further de¬ 
scribe the shape of the intrinsic magnetization density for 
an antiferromagnetic domain wall and discuss its physical 
significance as a twisting of the homogeneous spinless AFM 
into a state with a finite spin. The intrinsic magnetization 
adds up in two- and three-dimensional extended domain 
wall systems and can affect the dynamics of domain walls 
subject to external magnetic fields and spin-polarized cur¬ 
rents. We discuss how these consequences can go beyond 
that of the purely quantum topological effects^^ observed 
in one-dimensional spin chains. 

Studies of domains in AFMs and descriptions of the 
shape and properties of antiferromagnetic domain walls 
date back several decades.However, most of the exper¬ 
imental evidence of such domains were restricted to stud¬ 
ies of AFMs in which the collinearity of the sublattices is 
broken due to Dzyaloshinskii-Moriya (DM) anisotropy. In 
these studies, when the DM field or the external field van¬ 
ishes so does the equilibrium magnetization of the AFM. 
Consequently, the detection of domain walls and their dy¬ 
namics in compensated AFMs remains an experimental 
challenge. However, antiferromagnetic domain walls are 
known to exist and have been experimentally observed, e.g., 
in monolayers of antiferromagnetic Fe,^^ in the elemental 
AFM Cr,^^ and in the antiferromagnetic insulator NiO.^^ 
Antiferromagnetic domains and domain walls can also be 
tailored by manipulating the ferrimagnetic precursor layer 
before cooling the AFM below the Neel temperature.^^ Ob¬ 
servation of individual domains in AFMs can be done, e.g., 
using X-ray magnetic linear dichroism.^^’^^ 

A key aspect of detecting the dynamics of antiferromag¬ 
netic domain walls is whether such solitons of staggered 
magnetic order are associated with a spatially constricted 
magnetization density. Ref. 20 argued that such a mag¬ 
netization exists and that the earlier studies of antiferro¬ 
magnetic spin chains missed certain parity-breaking terms 
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in the transition from the discrete spin model to the con¬ 
tinuum approximation. The antiferromagnetic Heisenberg 
Hamiltonian has been mapped to the non-linear a model 
for the continuous staggered order parameter.How¬ 
ever, in Haldane’s seminal work on large-spin Heisenberg 
AFMs,^^ no apparent parity-breaking terms survived the 
transition to the continuum model. In Haldane^s map- 
the continuum field that is conjugate to the 
antiferromagnetic order parameter describes the dynamic 
magnetization only (see Sec. HD). Using a slightly dif¬ 
ferent parametrization of the antiferromagnetic order and 
the magnetization field, Ivanov et later demon¬ 

strated that the energy functional based on the one¬ 
dimensional antiferromagnetic Heisenberg model indeed 
contains a parity-breaking term in the continuum limit and 
that this term must be taken into account to describe the 
equilibrium magnetization of a domain wall. The parity¬ 
breaking term included in Refs. 21 and 22 is not equivalent 
to the well-known ’’topological 0 terms”which arise 
in effective a model Lagrangians for one-dimensional anti¬ 
ferromagnetic spin chains and are responsible for quantum 
effects such as Haldane’s conjecture. 

The recently increased interest in AFMs as active spin- 
tronics components has spawned a number of effective mod¬ 
els for antiferromagnetic dynamics.These recent 
models mostly adopt the non-linear a model without intro¬ 
ducing a Hamiltonian that includes parity-breaking terms 
that lead to the intrinsic magnetization of antiferromag¬ 
netic textures. The absence of parity-breaking terms in 
these models may be due to different definitions of the con¬ 
tinuum fields, or these terms may have been disregarded in 
the transition to the continuum limit due to specialized 
symmetry requirements, which only hold for homogeneous 
AFMs. Whether the intrinsic magnetization of extended 
two- and three-dimensional systems can lead to qualita¬ 
tively new physics for the dynamics of antiferromagnetic 
textures under the influence of external forces remains an 
open question that we seek to address in this article. 

The intrinsic magnetization of antiferromagnetic tex¬ 
tures is small. A domain wall in a one-dimensional antifer¬ 
romagnetic spin chain exhibits intrinsic magnetization that 
is in total no larger than the spin of one sublattice^^’^^. It is 
therefore unlikely that such a small magnetic moment can 
be directly detected in the near future. However, the pres¬ 
ence of the small spin of domain walls in one-dimensional 
spin chains manifests itself through quantum effects. 

In higher-dimensional extended systems, such as synthetic 
AFMs, the magnetization of a textured multilayer may be 
of appreciable size^^. Furthermore, in thin films or in bulk 
AFMs, which is the focus of our study, the intrinsic mag¬ 
netization of a transverse domain wall is additive in the 
perpendicular directions. The result is a macroscopic mag¬ 
netization that can be more easily excited and detected and 
that can influence the dynamics of AFMs beyond that of 
purely quantum effects. 

The paper is organized as follows. In the next section 
(ii), we take the continuum limit of the Heisenberg Hamil¬ 
tonian, describe the origin of the intrinsic magnetization, 
and discuss the consequences for the antiferromagnetic dy¬ 
namic equations. We also compare our model to Haldane’s 
alternative mapping of the continuum fields. This compar¬ 


ison demonstrates that the continuum fields in these two 
parametrization procedures have critically different physi¬ 
cal interpretations. In Sec. HI, we describe the magneti¬ 
zation profile of a domain wall and discuss generalizations 
to higher-dimensionsal systems. We show how the intrinsic 
magnetization leads to qualitatively new physics and that 
the domain wall can be moved by a spatially inhomoge¬ 
neous magnetic field that couples to the intrinsic magne¬ 
tization. In Sec. IV, we present numerical results for the 
motion and control of an antiferromagnetic domain wall 
and show that we can create potential wells for the domain 
wall with spatially constricted magnetic fields. In Sec. V, 
we discuss the experimental consequences of the intrinsic 
magnetization for extended systems in 2D and 3D. Sec. VI 
concludes the discussion. 


II. THEORY 

Our starting point is the Heisenberg Hamiltonian due to 
the exchange coupling between classical spin vectors on a 
lattice^^ 

= J ^ , (1) 

where the positive exchange energy, J > 0, describes an an- 
tiferromagnetic ground state, (a, /3) denotes a sum over all 
nearest neighbor lattice sites described by the two sublat¬ 
tices a and /3, where each spin at a has Nn nearest neigh¬ 
bors of type /3, and vice versa, a and /3 are D-dimensional 
vectors, where D is the dimensionality of the AFM. We 
proceed by describing the simplest model, the D = 1 anti¬ 
ferromagnetic linear spin chain with easy-axis anisotropy, 
and later generalize our results to 2D and 3D in Appendix 
A. The focus of our subsequent sections are on extended 
3D AFMs in which the order parameter varies along one 
dimension only. 

A. Free energy functional for ID 

We consider a linear spin chain with 2N atomic lattice 
sites, where the spins on half of the lattice sites, denoted 
by Of, minimize their energy by pointing in the opposite 
direction of the spins on their Nn = 2 nearest-neighbor 
lattice sites, denoted by /3, and vice versa. For the AFM, 
we impose the boundary conditions that the spin on the 
left end of the spin chain is of type a, whereas the right 
end of the chain is occupied by a /3 site. Therefore, in the 
ground state, the AFM is fully compensated, and the total 
spin vanishes. We define the 2 ; axis as the magnetic easy 
axis. The classical Heisenberg Hamiltonian including the 
easy-axis anisotropy is 

N,N / N N \ 

= J E ^ b (2) 

(a,/3> \ a g J 

where K is the anisotropy energy. In typical easy-axis 
AFMs, the exchange energy dominates, \J\ ^ \K\. The 
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classical ground state of the Hamiltonian (2) is degenerate, 
(Sc^, 8 / 3)0 ^ ±(5'^, —S'i), where S (in units of h) is the 
spin on a single atomic lattice site. 

We now introduce the standard definitions (see Sec. IID 
for a comparison with an alternative definition that is oc¬ 
casionally mistaken to be equivalent to the present model) 
of the magnetic and staggered order parameters, and 
li, on a two-sublattice linear lattice parametrized by i: 


h = 


SL + 

2S ’ 

’ 


(3a) 

(3b) 


where we have paired the sublattice spins and at 
unit cell i running over a total of N antiferromagnetic unit 
cells. In this convention, m? + 1 ? = 1 and the spins in unit 
cell i can be expressed as follows: 

Sj, = 5(m, + 1,), (4a) 

= 5(m, - 1,). (4b) 


After introducing the magnetization vector 111 ^ and the 
staggered order parameter 1^, the Heisenberg Hamiltonian 
( 2 ) reduces to a sum over antiferromagnetic lattice points: 


AT-l 

Hid = JS^ - li)[(mi + 1*) + (mj+i + Ij+i)] 


+J52(m2,-1^) 

N 

-KS"^ [(mi + \i)l + (mi - \i)l] . (5) 

i 


We continue by using the identities = m? + 

- (m^+i - m^)^ and (liOi^+i - mj^+i) = li(mi+i - 
m^) — mi ( 1^+1 — li) to rewrite the bulk part of Eq. (5) as 
follows: 

N 

H,d « 2JS^ - 1?) 

i 

T q2 

+ ^ E [(*»+! “ 

i 

AT-l 

+JS‘^ [mi(li+i - li) - li(mi+i - mi)] 

N 

+ (6) 


where we have disregarded the vanishingly small energy 
contribution — J5'^(m^ -|- — n |^)/2 from the unit 

cells at the edges. 

Next we go to the large-A" limit and take the con¬ 
tinuum approximation, allowing us to write Aid ~ 
/(di/A)'H(l, l',m,m' ), where A is the length of the antifer¬ 
romagnetic unit cell and 1' and m' are the (dimensionless) 
spatial derivatives of the staggered field and the magneti¬ 
zation, respectively, di is an infinitesimal length element 
along the spin chain. For A = 1 , A = 2d^ where d is the 


H 


a, |3 


H 
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FIG. 1. For a simple linear spin chain with antiferromagnetic 
exchange coupling, the Heisenberg Hamiltonian is not invariant 
under sublattice exchange, 8^ S^, if the order parameter 

is spatially inhomogeneous. (Top) A simplified sketch of a 6- 
spin 90° texture. Exchanging the spins on sublattices ol and /3 
(bottom) creates a more disordered phase that costs additional 
exchange energy, hence Ack ,/3 < In the continuum limit, 

this energy difference is captured by the parity-breaking term 
in the antiferromagnetic energy functional. 


nearest-neighbor spacing in the linear chain. The energy 
density (apart from a constant and in units of energy) is 

2210 ( 1 , l',m, m') = JS'^[4|mp + jl'p - jm'p 
+(m • 1' - 1 • m')] - KS"^ [(1 -zf+ (7) 

We note that the fourth exchange term in Eq. (7) has an 
unusual parity-breaking form^^’^^ and is an odd function 
of the order parameter 1 . 

In the models of AFMs that we consider, the two- 
sublattice linear lattice in ID, the centered squared lat¬ 
tice in 2D, and the body-centered cubic lattice in 3D, the 
Heisenberg Hamiltonian is not invariant under sublattice 
exchange {a f3) if the order parameter is spatially in¬ 
homogeneous, see Fig. 1. However, there is an ambiguity 
in the pairing of spins 8 ^^, and 8 ^ and the definition of the 
order parameter 1^ in Eq. (3b). One might as well choose 
1 - = — 1 - as the order parameter, and consequently, one usu¬ 
ally demands that the bulk Hamiltonian is invariant under 
the transformations 1 ^ ^ — 1 ^ and ^ because the 
two possible choices of the order parameter are physically 
equivalent. Under these transformations, the definitions of 
8 J^ and 8 ^ in Eqs. (4) also change, and the fourth exchange 
term in Eq. (7) undergoes an additional sign change. The 
energy functional is therefore invariant with respect to the 
two equivalent definitions of the order parameter but not 
invariant under sublattice exchange. In the latter case, 
the ordering of the spins changes, resulting in a larger ex¬ 
change energy penalty for inhomogeneous AFMs. A sim¬ 
plified sketch of this energy difference is shown in Fig. 1 
for a 6 -spin chain with a 90° texture. In the bottom spin 
chain the a and /3 sublattices have been exchanged, leading 
to a more disordered phase that costs additional exchange 
energy. This result generalizes to an arbitrary number of 
spins in a linear textured spin chain. 

To describe the order parameter dynamics of the AFM, 
it is useful to work in the exchange approximation^^, | J| ^ 

I A I, and consider slowly varying antiferromagnetic tex¬ 
tures. In this case, |mp <C |lp, and we can disregard terms 
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that are of higher order than |mp, such as the magnetic 
anisotropy energy term and the magnetic stiffness term in 
Eq. (7). We choose the spin chain axis to be along the 
^ axis and introduce the normalized staggered vector field 
n(z,t) = We can consequently write the en¬ 

ergy density as a function of the deviations dzU {= dnjdz) 
and m from the ground state. After integrating by parts, 
we arrive at the free energy density for the linear antiferro¬ 
magnetic spin chain to the lowest order in deviations from 
an equilibrium state: 

'HiD(n,92n,m) = ^\m\'^ + ^\dzn\'^+L{m-dzn)-^{n-zf . 

(8) 

The equation has the following parameters: the homoge¬ 
neous exchange energy a = the exchange stiffness 

terms A = /S?JS^ and L = 2AJ5'^, and the anisotropy 
energy Kz = 2KS‘^. Here, a finite L lifts the degeneracy of 
the sublattice exchange. 


B. Free energy functional for D > 1 


In Appendix A, we generalize the free energy of Eq. ( 8 ) to 
2D and 3D for the centered squared and the body centered 
cubic unit cell, respectively. We find that the generalized 
free energy density in the exchange approximation is given 

by 


7^(n, m) 



i i^j 


+LY^{ni-din)-^{n-zf, (9) 


where a = 4 NdJS^, A = Nd^^JS^/2, L = NdAJS^, 
Kz = 2KS^^ and Nd is the number of nearest neigh¬ 
bors. Ni = 2, ^2 = 4 for the squared lattice, and 
depends on the choice of unit cell, 6 for the simple cubic cell 
and 8 for the body-centered cubic cell. The stiffness part 
of the above Hamiltonian density contains two apparent 
anisotropic terms, ^ ( 9^11 • dju) and ^ (m • diu). However, 
in the following, we show that after eliminating the degrees 
of freedom associated with m, the effective Lagrangian re¬ 
duces to the non-linear a model and the resulting antifer¬ 
romagnetic spin-wave dispersion remains isotropic . 

This isotropic dispersion is in contrast to the anisotropic 
dispersion relation resulting from the exchange term iden¬ 
tified by Lifshitz and Pitaevskii^^, which is similar but 
not identical to the third term in Eq. (9). Lifshitz and 
Pitaevskii consider only the small deviation (n ^ 
no + n^) from the equilibrium homogeneous antiferro¬ 
magnetic spin configuration and add the exchange term 
^ (m • dzn± — nx • dzm) to the free energy density. Com¬ 
pared to Eq. (7), this also results in a surface anisotropy 
^ (no • dzin), which (after integration over the space) fa¬ 
vors magnetization build-up on the edges of the AEM. 
Consequently, the dispersion relation for this model is 
anisotropic. The parity-breaking exchange term (^ L) in 
the above free energy density (7) differs from the term of 
Lifshitz and Pitaevskii because it involves n rather than 
and does not violate the isotropic dispersion relation 


of antiferromagnetic spin waves due to small variations in 
the staggered field n. This is also the case for D > 1. 
Neglecting the parity-breaking term as being of leading or¬ 
der in the exchange energy would imply that an AEM at 
equilibrium exhibits no intrinsic magnetization, even when 
textures in the staggered field are present. 


C. Lagrangian density and equations of motion 


The equations of motion for the staggered field n and 
the magnetization field m can be found from, e.g., linear 
combinations of the equations of motion for the sublat¬ 
tice spins Sq, and Equivalently, we may proceed by 

constructing the Lagrangian density and directly compute 
the dynamic equations for n and m from the variation of 
the Lagrangian with respect to these fields. Our starting 
point is the generalized free energy density in the exchange 
approximation, Eq. (9). The Lagrangian density can be 
constructed as £ = /C — H, where 1C is the kinetic energy 
term. Analogous to the procedure for constructing the ki¬ 
netic term for a single spin in a ferromagnet,^^’^^ 1C can be 
constructed from the Berry phase of the spin pair Sq, + 
that constitutes the antiferromagnetic unit cell: 


J icdv = -Sh 


^ ^ ’ Sq + ^ ^ , 

oc /3 


( 10 ) 


where it is convenient to choose the gauge potential Aq^^) 
such that the spin-pair Berry phase vanishes in the strictly 
antiparallel configuration, Sq = — S/ 3 . One such choice is 
Aq(/ 3 ) = — 0 ck(/ 3 ) cos 6 >q(^)/ sin 6 >Q(^) in the spherical coor¬ 
dinate system, where 0 is the polar angle and 0 is a unit 
vector along the azimuth. This gauge is identical to that 
which is normally used to describe the kinetic energy of 
a single spin in ferromagnets^^ and generalized to a two- 
sublattice model with antiparallel spin configuration. By 
expanding the spin pair Berry phase in small deviations 
from the antiparallel configuration, 6^/3 tt — ( 6 >q + 60) and 
0/3 ^ TT + ( 0 Q + (i 0 ), and transferring back to the [n, m] 
basis, the kinetic term in the continuum approximation is 
given by^^’^^ 


1C = pm(h X n), (11) 

where p = 28% is the magnitude of the staggered spin 
angular momentum per unit cell and we have disregarded 
terms of the order |m 0 and higher. 

Varying the Lagrangian with respect to the magnetiza¬ 
tion m and the staggered field n gives the coupled Landau- 
Lifshitz equations of motion 

h = cJm X n , ( 12 a) 

ih = cJn X n + cJm X m , (12b) 

where damping is typically phenomenologically 
introduced.In the transverse basis, where |np = 1, no 
term of the form ^ (cJn x m) (as present in the dynamics 
of 1 in, e.g.. Ref. 10) appears in Eq. (12a), which is valid 
in the exchange approximation and includes terms up to 
second order in small deviations from equilibrium. The 
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effective magnetic and staggered fields (in units of s are 
defined as functional derivatives of the total free energy 
H: 

pi^m = -|^ = -am - , (13a) 

= A(V^n + 

Kz{y\ ' z)z ^ (13b) 

where we have defined the sum over spatial derivatives 
in all directions as d'^ = — 

9^/ (didj). In the Appendix (A) we discuss how these 
anisotropic differential operators arise in 2D and 3D. 

In the absence of external forces in the effective magnetic 
field, Eqs. (12a) and (13a) give^^ 

m = -h X n - —d^n , (14) 

a a 

which indicates that the magnetization field m is simply 
a slave variable that follows the temporal and spatial evo¬ 
lution of the staggered field n. We note that if we ne¬ 
glect the parity-breaking term in the free energy (L ^ 0), 
the intrinsic magnetization of a textured AFM vanishes 
at equilibrium. Our analysis shows that for our partic¬ 
ular parametrization of the continuum fields, this parity¬ 
breaking term is an important part of the transition from 
the discrete spin model to the continuum approximation 
and cannot be disregarded. 

Equation (14) allows us to eliminate m and write an 
effective Lagrangian density for the staggered field and its 
derivatives as 

£(n,h,ain) = ^\h\^ ^ E 

i 

+ n) + 

a A 

I 

This Lagrangian density describes the anisotropic non¬ 
linear cr model with a kinetic topological term (third 
term).^^’^^’^^“^^ This topological term is a by-product of 
the elimination of m from the Lagrangian. It can be shown 
that this term has the form of a total derivative^^. Conse¬ 
quently, it has no effect on the effective equations of motion 
for n or the domain wall dynamics that we describe in the 
next sections. We will not discuss in any detail the quan¬ 
tum effects of the topological term in the following. 

D. Comparison with Haldane’s mapping 

We digress for a moment to compare the one-dimensional 
model described above with a commonly used alterna¬ 
tive definition of the continuum fields known as Haldane’s 
mapping^^’^^’^^ of the antiferromagnetic order parame¬ 
ter. We include this comparison because the different 
parametrizations are not equivalent and are, therefore, re¬ 
current sources for confusion. In contrast to the Hamil¬ 
tonian approach described by Eqs. (3) and (4), Haldane’s 
parametrization maps each spin in the spin chain at cite i 
onto two continuum fields: 

Sj/S' = - ihf + ihj, (16) 



[iij, mj 

i 

1 


/ \ 

t \ 

Hamiltonian approach: 


Ci+1 Ci+1 

^ a ^ |3 

- t 4 t T. 

Haldane’s mapping: 

S2i.l 

^21+1 ^24+2 


FIG. 2. In the Hamiltonian approach (top), Eqs. (3) define 
values for the staggered field and the magnetization field m^ 
at the center of every antiferromagnetic unit cell labelled by i. 
In Haldane’s mapping (bottom), every single spin is mapped 
onto two continuum fields, the Neel field h and the ’’canting” 
field m. 

where h is the unitary Neel field and m is the ’’canting” 
field. We note that this mapping introduces extra degrees 
of freedom, which must subsequently be reduced by limit¬ 
ing the Eourier components of the fields n and m to include 
only long-wavelength excitations.^^ 

Figure 2 compares the labelling of the spins in the Hamil¬ 
tonian approach used in this work with that of Haldane’s 
mapping. By equating the expressions for and in 
Eqs. (4) and their corresponding expressions in Haldane’s 
parametrization, we find the relationship between the con¬ 
tinuum fields in the two different parametrizations: 

nij + Ta-i\Jl - mf = + m2i_{17a) 

rrij - Ta.i\Jl - mf = - mij + m2i (17b) 

In the exchange approximation, m <C n and m <C h. Keep¬ 
ing only the lowest order contributions in the magnetization 
m and the canting field m, it follows that 

“i ~ -^(n 2 i-i + n 2 i) + ^(m 2 i-i - m 2 i), (18a) 

“li ~ -^(n2i-i - n2i) + ^(m2i-i + m2i), (18b) 

where we have disregarded terms of the order |mp and 
|ihp and higher. 

For small-angle spatial variations in the continuum fields, 
we use the gradient approximation to find the field val¬ 
ues for n and m at the center of each unit cell: h^^i /2 ~ 
fii + (A/4)5;2hi and m^+i /2 ~ ixii + {/S./A)dzmi, where 
A/2 = d is the nearest neighbor distance and n(m)^+i /2 
represents the Neel (canting) field at the midpoint between 
the spins and S^+i. Inserting these lowest order gradi¬ 
ent approximations into Eqs. (18) results in a one-to-one 
relationship between the continuum fields of the Hamil¬ 
tonian approach and Haldane’s parametrization. Corre¬ 
spondingly, the mapping between the two different repre¬ 
sentations reduces to n ^ —n + {/S./A)dzm + 0(|mp) and 
m ^ m — {/S./A)dz^ + 0(|mp). 

It is critically important that the continuum fields h and 
m of Haldane’s mapping are not identical to the staggered 
and magnetization fields n and m used in the present work. 
By inserting the mapping between the two parametriza¬ 
tions into the energy functional in Eq. (8) and keeping only 
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terms of the order |mp in the exchange approximation, we 
find the continuum limit energy functional of Haldane’s 
mapping: 

^Hai(n,5^n,m) = ||mp + ^|5^np- . (19) 

This result conclusively shows that the parity-breaking ex¬ 
change term in Eq. (8), which is a result of the proce¬ 
dure of breaking the lattice into spin pairs, vanishes af¬ 
ter a suitable transformation of the continuum fields, e.g., 
m ^ m — (A/4)9;2h. In other words, when applying Hal¬ 
dane’s mapping procedure, the parity-breaking exchange 
term does not appear in the energy functional. An overall 
requirement, however, is that the physics remains the same, 
including the existence of the intrinsic magnetization. 

Although the Hamiltonian approach used in the present 
work and Haldane’s mapping are both valid continuum 
representations of spin systems with antiferromagnetic ex¬ 
change coupling, a crucial difference exists for the physi¬ 
cal interpretations of the continuum fields, which are not 
equivalent in the two representations. The equilibrium 
value of the canting field m of Haldane’s mapping vanishes, 
also when h is inhomogeneous. Therefore, m represents the 
dynamic magnetization induced by temporal variations of 
the order parameter h and not the total magnetization. 
Consequently, the coupled equations of motion for h and 
ih are not of the same form as Eqs. (12) and (13). In par¬ 
ticular, the expression for the canting field: m ^ h x h, 
which is analogous to Eq. (14), does not include a term 
proportional to the gradient of h. This fact may be an 
important reason why the intrinsic magnetization is easily 
missed in models based on Haldane’s mapping. 

In the Hamiltonian approach, on the other hand, m can 
be interpreted as a magnetization density in the sense that 
the total accumulated spin (both intrinsic and dynami¬ 
cal) of the AEM can be found from integration, M/S' = 
J m dV. Eor antiferromagnetic textures, this integral is 
generally nonzero even for static spin systems. Although 
the canting field m in Haldane’s mapping does not in¬ 
clude the intrinsic contribution to the magnetization den¬ 
sity, the total spin can instead be found from the relation 
M/S ^ The intrinsic magneti¬ 

zation can be identified as arising from the first terms in 
the sum. Eor a slowly varying h in, e.g., the z direction, 
~ [nzizi)-nz{z 2 N)]/ 2 y^ which is gener- 
ally nonzero for textured AEMs. In the following analysis, 
we continue using the Hamiltonian approach, in which the 
continuum field m is interpreted as the total magnetization 
density. 

E. Antiferromagnetic spin waves and spin current 

To study small harmonic excitations from a homogeneous 
AEM, we construct the effective equation of motion for the 
staggered vector field n by combining Eqs. (12a) and (12b) 
while retaining the constraint |np = 1: 

n X (h X n) = x [{aA — I/^)V^n aKz{n • z)z] x n. 

^ ( 20 ) 


The parity-breaking exchange term leads to the renormal¬ 
ization of the exchange stiffness A ^ A* = (A — L^/a) = 
A/2 but otherwise leaves the equation of motion (20) invari¬ 
ant in linear response.The topological term in Eq. (15) 
has no effect on the effective equations of motion for n, as 
expected. 

Insertion of a small harmonic excitation from the ground 
state in time and space, n(r, t) z-\-Sn± exp [i{ujt — k • r)] 
into Eq. (20) results in the usual ’’relativistic” antiferro¬ 
magnetic dispersion relation 

= ^[aA*k^ + aK,], (21) 

P 

where k = |k|. In the isotropic limit, Kz 0, which results 
in the familiar linear dispersion 

uoi = ck , (22) 

where c = Nj^SJA/{2h) is the spin wave phase veloc¬ 
ity. Eor A = 2d/^/D^ where d is the nearest-neighbor 
distance, and for hypercubic lattices, where Nn = 2D^ 
Eq. (22) agrees with Eqs. (13) and (20) in the semi-classical 
treatment in Ref. 40, as well as with Holstein-Primakoff 
calculations^^’^^ and Haldane’s D = 1 result.We note 
that the parity-breaking term (^ L) does not lead to an 
anisotropic dispersion relation, such as the term in Lifshitz 
and Pitaevskii.^^ On the contrary, the inclusion of such a 
term is important to arrive at the correct dispersion rela¬ 
tion in the classical continuum limit. 

The intrinsic magnetic moment of antiferromagnetic tex¬ 
tures will necessarily influence how spin currents in inho¬ 
mogeneous AEMs are described. A continuity equation for 
the spin angular momentum transfer in the AEM caused by 
the exchange interaction can be constructed from Eq. (12b) 
as pm -j- ^ • di^s,i = 0- The spin current polarized along i 
is 

Jg i = A*5in X n — —h, (23) 

’ a 

where we have used Eq. (14) to eliminate m. Equation (23) 
explicitly shows that a time-varying antiferromagnetic tex¬ 
ture is equivalent to spin angular momentum transfer, a 
relationship that easily can be missed by models for the 
staggered dynamics that disregard the intrinsic magneti¬ 
zation. This result may have implications for antiferro¬ 
magnetic spin pumping from textures^^ because the col¬ 
lective motion of the antiferromagnetic order parameter is 
equivalent to a current of spin angular momentum. In one¬ 
dimensional textures, pL/a = Shd, thus indicating that 
textures that oscillate at frequency T~^ produce a spin- 
current corresponding approximately to a single spin mov¬ 
ing one lattice spacing per period of oscillation T. 

F. Consequences for staggered dynamics 

In effective models for the dynamics of the staggered vec¬ 
tor field n, the magnetization field m plays the role of a 
slave variable that follows the temporal and spatial evolu¬ 
tion of n. When no external forces couple directly to the 
intrinsic spin in the AEM, the parity-breaking term in the 
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energy functional L) only leads to a renormalization of 
the exchange stiffness and has no other effect on the dy¬ 
namic equations. However, we show in the following that 
by including external magnetic fields or spin-polarized cur¬ 
rents, the dynamics of the antiferromagnetic order param¬ 
eter can also be altered indirectly through the excitation of 
the magnetization density field m. 

The spin-transfer torque on ferromagnetic textures is 
normally considered a second-order effect in ATMs when 
acting only on the small magnetization m(t) induced by the 
time variation of the staggered field, h. If ATMs also ex¬ 
hibit intrinsic magnetization, the spin-transfer torque from 
spin-polarized currents on the magnetization m may be¬ 
come more important. However, because the magnetiza¬ 
tion is first order in the spatial variation of the staggered 
field, m ^ the Berger spin transfer torques (Eqs. (5) 
and (6) in Ref. 52) are of the order J smaller than the 
driving forces acting directly on textures in the staggered 
field, first identified in Ref. 16. In this case, the intrinsic 
magnetization of ATMs will lead to higher-order correc¬ 
tions to the current-induced torques that couple directly 
to the staggered field. In antiferromagnetic thin films with 
strong surface anisotropy or in special cases in which a 
strained geometry suppresses the torques on the staggered 
field, the Berger torques on the textured magnetization 
could become important. We will not discuss the effects 
of spin-polarized currents any further in the following. 

Instead, we focus on the effect of an external magnetic 
field H that couples directly to the intrinsic magnetiza¬ 
tion of antiferromagnetic textures. To illustrate this phe¬ 
nomenon, we add the Zeeman interaction to the free energy 
density, Hr = H — P 7 (H • m), where 7 is the gyromagnetic 
ratio. The external magnetic field induces a small magnetic 
moment density in the AFM, and the magnetization field 
m is altered according to 

m = -h X n — —d^n + x (H x n), (24) 

a a a 

where the cross products enforce the constraint n • m = 0 . 
Inserting this result in the Lagrangian gives the effective 
Lagrangian density for an AFM under the influence of an 
external magnetic field H: 

^ la*-" “ ^ ~ T 

i 

+ ^(n ■ z)^ + {n X h) 

A a 

I 

(25) 

i 

This Lagrangian density agrees with that proposed in 
Ref. 49, with the exception of the second to last topolog¬ 
ical term and the last term, which couples the external 
magnetic field and textures in the antiferromagnetic order. 
In the following, we show how this coupling between mag¬ 
netic fields and the gradient of the staggered field allows 
the movement of domain walls in ATMs to be controlled 
by spatially varying magnetic fields. This result has not 
been reported previously. 

Utilizing the method of collective coordinates^^’^^, we as¬ 
sume that the temporal dependence of the staggered vector 


field n(r, t) is held by a set of coordinates {aj{t)} that de¬ 
scribe the time evolution of textures in the AFM, such that 
n(r, {aj{t)}). In this case, the time derivative of the stag¬ 
gered field can be written as h = djdajn. We earlier 
demonstrated that in AFMs, the collective coordinates can 
be viewed as quasi-particles with an effective mass reacting 
to external forces and following Newton’s 2nd law.^^ The 
equation of motion for the collective mode aj is 

• • ncy 

= (26) 

where /a) f dV(9ain • Sajii) is the effective mass, 

a is the phenomenological Gilbert damping parameter for 
AFMs, and F* are the forces that act on the collective 
excitations. F* = + F^^^ can be split into the internal 

exchange forces F^^^ = da^H^ which are derivatives of the 
free energy with respect to the collective modes, and the 
external forces F^^^. We focus here on an external magnetic 
field as the only external force that acts on the AFM, giving 

= ^J dV[pil ■ (n X da,n) + L{da,n ■ 5*H)], (27) 

where, in addition to the previously identified reactive force 
on the collective coordinates in AFMs due to time-varying 
magnetic fields,we now identify a new force induced by 
a spatially inhomogeneous magnetic field. This force will 
necessarily influence how antiferromagnetic textures are ex¬ 
cited by external magnetic fields. 

III. DOMAIN WALL DYNAMICS 

In this section, we return to systems where the order 
parameter varies along one dimension and discuss how the 
intrinsic magnetization influences the motion and detection 
of solitons in quasi-one-dimensional AFMs. Although the 
texture is assumed to vary only along one direction, the 
nearest neighbours to each spin may also exist along two 
(2D) or three (3D) axes. Later, we show how a Neel domain 
wall can be accelerated and controlled by a stationary and 
spatially inhomogeneous magnetic field. 

A. Antiferromagnetic domain walls 

In one-dimensional spin chains, the spatial variation of 
the staggered field n is constricted to the spin chain axis, 
z. At equilibrium, the time evolution of the staggered field 
and the magnetization vanishes, h = 0 and ih = 0, and 
Eq. (20) gives 

no X [A*9^no + Kz{no • z)z] x no = 0 . (28) 

By introducing spherical coordinates for the 
normalized staggered vector field as no ( 2 ;) = 
[sin Oq cos 00 5 sin Oq sin 0o 5 cos Oq], a series of solutions 
for the above equation can be found from 

9;s0o = O, (29a) 

d^Oo = ^ sin Oo cos Oq , (29b) 




FIG. 3. Sketch of the intrinsic magnetization m(z) (red) (not 
to scale) of one-dimensional (a) Neel and (b) Bloch (not used 
in the calculations) domain walls in the order parameter 11 ( 2 :) 
(green). The equilibrium magnetization profile was calculated 
from Eq. (14). We note that the magnetization is so small 
that in a one-dimensional system, the domain wall spin must 
be treated quantum mechanically. However, for higher dimen¬ 
sional extended systems, the total spin of domain walls could 
be of appreciable size because the intrinsic magnetization is ad¬ 
ditive in the perpendicular directions. 


where A = The trivial solution to Eqs. (29) 

is Oo{z,t) 0, which corresponds to a homogeneous 
AFM where all the spins are polarized along the posi¬ 
tive/negative z direction. The excited state is given by 
Oq = 2 arctan[exp(z/A)], the Walker domain wall.^^ In 
this Neel configuration, sin6>o = ±sech(z/A) and cos6>o = 
±tanh( 2 :/A), which ensures that iiq = 1. A is the half¬ 
width of the domain wall. Inserting the results from the 
Heisenberg model, we find that the domain wall half-width 
A = is given by a competition between the ex¬ 

change and anisotropy energy scales, as expected. 

The intrinsic magnetization associated with the antifer¬ 
romagnetic domain wall at equilibrium is given by Eq. (14) 
when n = 0: 

sech( 2 :/A) tanh( 2 :/A) cos 0o 

sech(2:/A) tanh(^/A) sin^o , 

—sech^(2;/A) 

(30) 

where the sign determines whether the Neel domain wall is 
head-to-head or tail-to-tail. The magnetization profile of a 
head-to-head Neel domain wall and the profile of an anti¬ 
ferromagnetic Bloch domain wall are presented in Fig. 3. 
The total magnetic moment in the z direction contained in 
a head-to-head domain wall configuration is^^’^^ 

^ J dz{mo • i) = S', (31) 

This result demonstrates that domain walls in the anti¬ 
ferromagnetic order induce a finite magnetization propor¬ 
tional to the spatial derivative of the staggered field and 


mo =-S^no = ± — 

a 2A 


that the direction of the magnetization depends crucially 
on the boundary conditions of the AFM, e.g., in the case 
of the Neel wall whether it is head-to-head or tail-to-tail. 
This result is intuitively easy to appreciate: because both 
edge spins (at an a and P site) point in the same direction, 
the 180° twist turns the homogeneous spinless AFM into 
a spin-S' object. The domain wall is a non-linear excita¬ 
tion of the homogeneous AFM and carries the spin S. The 
creation of a domain wall can therefore be interpreted as 
a twisting of the homogeneous spinless AFM into a con¬ 
figuration with a finite spin S that is located around the 
domain wall center. 

A consequence of the intrinsic magnetization of domain 
walls in one-dimensional spin chains is that for AFMs with 
half-integer S, the ground state, which is doubly degen¬ 
erate, occurs for stationary domain walls,and not for 
processing domain walls, as predicted in Ref. 26. Another 
consequence is that the motion of domain walls in AFMs 
is equivalent to spin angular momentum transfer, as con¬ 
firmed by Eq. (23). The identification of antiferromagnetic 
domain walls as single-spin carriers may become important 
for future applications in antiferromagnetic spintronics. 


B. Domain wall motion 


We consider a (slowly) moving tail-to-tail domain wall 
profile n[^,a^(t)] corresponding to the dynamic soliton so¬ 
lution 0{z,t) 2arctan{exp[(^ — (t))/A]} and (p{t) -)> 

The domain wall shape is assumed to be rigid, so 
that the temporal dynamics are held by the collective coor¬ 
dinates {an{t)} the domain wall tilt an¬ 

gle with respect to the x-z plane and the position of the do¬ 
main wall center, respectively. Dissipation in AFMs is typ¬ 
ically added in a phenomenological mannerand is nat¬ 
urally incorporated in the collective coordinate approach. 
We add to the system a spatially varying magnetic field in 
the z direction, H = {0,0^ Hz{z)}. To the lowest order 
in the small external field and the velocities and 
we find that pw vanishes (although a constant precession 
(j)yj ^ ^ allowed in one-dimensional easy-axis systems) 
and that the domain wall center coordinate is accelerated 
according to 


aa . 
T 

P 


_7T 

7rpA 


H: 


int 


(32) 


where a is the dimensionless Gilbert damping parameter 
of the AFM. Depending on the spatial profile of the mag¬ 
netic field in the vicinity of the domain wall, the center 
coordinate will feel a force. The integrated magnetic field 
contribution is 



sech( —-— )tanh( —-— )Hz [z) , 

A A 


(33) 


where any non-even profile Hz{z) around the domain wall 
center coordinate gives rise to a finite acceleration of 
the domain wall. A homogeneous magnetic field does not 
accelerate the domain wall. In the steady state, the domain 
wall velocity saturates at / {'KaaX). We note 

that the domain wall velocity depends on the spatial dis¬ 
tribution of the external magnetic field. This dependence 
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opens up the possibility that nanoscale magnetic probes 
can accurately control the position of domain walls in, e.g., 
antiferromagnetic nanowires. In particular, a spatially con¬ 
stricted magnetic field can act as a potential well for the 
domain wall. In two-dimensional antiferromagnetic thin 
films, a spatially concentrated magnetic probe may attract 
spins from the edges of the AFM to form vortex states, see 
Sec. V. 


IV. NUMERICAL RESULTS 

To conceptually test the effect of a spatially inhomo¬ 
geneous magnetic field on the dynamics of an antiferro¬ 
magnetic domain wall, we have conducted numerical sim¬ 
ulations of generalized versions of Eqs. (12a) and (12b) in 
which we have phenomenologically included dissipation as 
in Ref. 13. We write the equations of motion in dimen¬ 
sionless form by scaling the time axis by t = pj and the 
spatial axis by the nearest-neighbor distance d. We solve 
the dimensionless equations of motion using the numerical 
method of lines with an adaptive time control. The system 
size is 2 : e [—500, 500] with the boundary conditions that 
n^(-500) = -1 and n^(500) = 1. 

Although domain walls in insulating AFMs, such as 
NiO, are approximately 100 nm wide,^^ we consider here 
the much shorter and more technologically important do¬ 
main walls observed in antiferromagnetic Fe-monolayers on 
W(OOl)^^, for which the geometric anisotropy is consider¬ 
ably larger. In such systems, the domain wall widths are 
only a few lattice spacings and the intrinsic magnetization 
is therefore relatively more important. For spin-1/2 parti¬ 
cles, for which the anisotropy energy per atom is 2.4 meV,^^ 
the time unit t 1 ps, the velocity unit v = dji ^ 300 ms“^ 
and the external field unit h = p^/Kz ~ 0.3 T. 

Fig. 4 presents the motion of a domain wall with half¬ 
width \ = Ad due to a constant magnetic field gradient. 
Because the domain wall spin in this particular Neel do¬ 
main wall is —S', the wall drifts toward lower magnetic 
fields to minimize its energy. The domain wall quickly 
reaches a steady-state velocity of approximately 50 ms“^. 
Fig. 5 presents how spatially concentrated magnetic fields 
can control and pin the position of the domain wall. By 
switching the pinning potential from the left to the right 
side of the domain wall, the position of the wall can be ac¬ 
curately controlled. The velocity of the center coordinate 
reaches more than 100 ms“^, and the transition from the 
left to the right pinning potential occurs in less than 100 
ps. 



FIG. 4. (a) Magnetic field strength as a function of distance 
along the spin chain. The magnetic field has a constant gra¬ 
dient of approximately 0.4 mT per lattice constant, (b) Time 
evolution of an antiferromagnetic domain wall moved by the 
magnetic field gradient. For clarity only the region z e[—25, 25] 
is shown. The domain wall slows down due to the finite dissi¬ 
pation when it reaches the region of the homogeneous external 
field. The maximum magnetic field strength is 10 mT, 

and the Gilbert damping constant is set to a = 0.01. 



V. HIGHER DIMENSIONAL EXTENDED 
SYSTEMS 

In this section, we discuss the possible experimental con¬ 
sequences for higher-dimensional textured systems, which 
typically extend in one or two perpendicular directions to 
the texture gradient axis. In such systems, the intrinsic 
magnetization can add up to a macroscopic number that 
is much larger than the spin on one atomic site. We also 
discuss the intrinsic magnetization of vortex states of the 


FIG. 5. (a) Magnetic field potential wells as a function of 

distance along the spin chain, (b) Time evolution of an an¬ 
tiferromagnetic domain wall controlled by the magnetic field 
potential wells. At t = 0, the domain wall is attracted to¬ 
ward a potential well at = —10 due to a spatially concen¬ 
trated magnetic field in the £ direction with the spatial profile 
Hz = Hosech[{z — Z-) / 10]z. In the time interval t = 140 ^ 160, 
the potential well to the left is turned off, and a similar mag¬ 
netic field-induced potential well is turned on to the right at 
= 10. 
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staggered order, which are two-dimensional analogs of the 
domain wall in the one-dimensional spin chain. At the end 
we briefly discuss the effects of pinning sites on the domain 
wall dynamics. 


A. Antiferromagnetic vortex states 

For D = 2 and in quasi-two-dimensional systems, such as 
antiferromagnetic thin films, non-trivial topological objects 
such as vortices^^ can form due to DM fields or external 
pinning. Fig. 6 shows the intrinsic magnetization m(x, y) 
associated with the spatially inhomogeneous staggered vec¬ 
tor field of such a two-dimensional object. The magnetiza¬ 
tion profile is calculated from Eq. (14). We note that the 
intrinsic spin of the vortex structure can be interpreted as 
a twisting of the spins in the homogeneous spinless AFM 
induced by spin rotations on the corners into a state with 
a finite spin located around the vortex core. The staggered 
vector field n(x,^) of this type of vortex structure is ro- 
tationally invariant around the vortex core along an axis 
normal to the x-y plane. The underlying spin structure, 
however, is not rotationally invariant, which is captured by 
the finite magnetization density m(x, y) of the vortex. The 
total spin of the vortex is S', as in the case of a domain wall, 
and the direction of the intrinsic spin depends crucially on 
the boundary conditions of the AFM, e.g., induced via ex¬ 
change bias pinning to ferromagnetic neighbors. 

The topological term in the effective Lagrangian density 
(15) for the staggered vector field n can possibly indirectly 
influence the dynamics of two-dimensional objects in the 
order parameter such as vortices or skyrmions. However, 
the complex two-dimensional dynamics of such topological 
objects are beyond the scope of this article and will not be 
discussed further. 


B. Extended domain walls in 2D and 3D 

Because the intrinsic spin of one-dimensional domain 
walls and two-dimensional vortices totals no more than the 
spin on a single atomic lattice site S', it is unlikely that the 
intrinsic magnetization associated with these antiferromag¬ 
netic textures can be reliably detected in the near future. 
Furthermore, to predict the correct excitation scheme of an- 
tiferromagnetic solitons, the intrinsic spin must be treated 
quantum mechanically because quantum fluctuations be¬ 
come important.^^ In higher-dimensional systems such as 
thin films or bulk AFMs, however, domain walls are not 
purely one-dimensional objects. Although the order pa¬ 
rameter can be defined as varying along one axis only, the 
nearest neighbors of each spin can exist along two (2D) or 
three (3D) perpendicular axes. In such systems, the in¬ 
trinsic magnetization of the domain wall accumulates over 
the total number of spin chains that constitute the domain 
wall structure. An example of the intrinsic magnetization 
of such an extended domain wall structure in, e.g., a nanos¬ 
trip is presented in Fig. 7. 

In bulk AFMs with domain structures in the order pa¬ 
rameter, the intrinsic magnetization forms planes along the 
domain boundaries. The total spin of these magnetization 


n(x,y) m(x,y) 
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FIG. 6. a) Sketch of a 2-dimensional antiferromagnetic vor¬ 
tex structure in the staggered order parameter n(x,y) (blue 
vector field). The portions of the staggered field pointing in 
the perpendicular direction close to the vortex core have been 
omitted for clarity, b) The resulting magnetization density pro¬ 
file m{x^y) (red vector field, not to scale) of the vortex state, 
calculated from Eq. (14). This intrinsic magnetization can be 
interpreted as a rearrangement of the spins on the corners so 
that the center of the vortex structure acquires a finite spin, c) 
A simplified sketch of a vortex structure with six spins along 
each edge ordered in centered squared unit cells with a (blue 
arrows) and 13 (red arrows) sites. Although the vortex struc¬ 
ture in the continuous staggered field appears invariant under 
the rotation of an arbitrary angle around the vortex core, the 
underlying spin structure is only invariant under axis inversion, 
[x, y] [—X, —y]. The total integrated spin of the vortex struc¬ 
ture is S, such as for the one-dimensional domain wall. The 
direction of this intrinsic spin is determined by the boundary 
conditions of the AFM. 


planes can be of appreciable size. In addition, for synthetic 
antiferromagnetic superlattices, in which the magnetiza¬ 
tion of each single ferromagnetic layer is much larger than 
S', the intrinsic magnetization associated with magnetic 
textures is accordingly larger and may be detectable. 


C. Effect of pinning sites on domain wall dynamics 

Pinning sites for domain walls can arise from impurities 
or crystal defects in the underlying lattice of AFMs. Al¬ 
though several studies have found that pinning effects in 
AFMs are small,dislocations and impurities diffuse 
the effects on a single spin. In quasi-one-dimensional spin 
chains the introduction of a single impurity atom can be 
enough to destroy long-ranged antiferromagnetic order and 
domain wall configurations. From such a perspective, the 
scenario studied in Sec. IV requires a perfect spin chain in 
strictly one-dimensional systems. However, because three- 
dimensional domain boundaries are typically sums of many 
one-dimensional spin chains, we expect the effects of pin¬ 
ning from impurities to be significantly smaller for domain 
wall systems that extend in the perpendicular directions 
than for one-dimensional spin chains. 
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FIG. 7. a) Sketch of a domain wall in the staggered order pa¬ 
rameter n(x, y) (blue vector field) in a two-dimensional AFM, 
e.g., a nanostrip. The domain wall configuration is repeated in 
the perpendicular direction, b) The magnetization vector field 
m(z,y) (red vector field, not to scale) calculated from Eq. (14). 
Each spin chain along the horizontal direction will contribute 
the spin Sx to the total spin of the two-dimensional domain 
wall structure. 


VI. CONCLUSION 


Starting from the discrete Heisenberg Hamiltonian 
with antiferromagnetic exchange coupling and easy-axis 
anisotropy, we have rederived the continuum limit of the 
free energy functional in the exchange approximation and 
conclusively shown that textures in the antiferromagnetic 
order exhibit intrinsic magnetization. In recent effective 
models for the dynamics of the antiferromagnetic order 
parameter, this intrinsic magnetization has been mostly 
disregarded. By comparing the Hamiltonian approach 
that we apply in this article with a commonly used al¬ 
ternative parametrization procedure called Haldane’s map¬ 
ping, we have shown that the continuum fields of the two 
parametrization procedures have crucially different phys¬ 
ical interpretations. As a result, the intrinsic magnetiza¬ 
tion can be easily missed in continuum models based on 
Haldane’s mapping. 

We have demonstrated that parity-breaking terms in the 
energy functional influence the dynamics of textured AFMs 
affected by external forces that couple directly to the intrin¬ 
sic magnetization. For extended domain walls in 2D/3D 
the influence of the intrinsic magnetization on the texture 
dynamics goes beyond that of the quantum effects observed 
for one-dimensional spin chains. By utilizing the method of 
collective coordinates, we have shown that a spatially inho¬ 
mogeneous magnetic field represents a reactive force on an¬ 
tiferromagnetic textures and can move a domain wall in an 
antiferromagnetic nanowire. This effect is directly linked 
to the intrinsic magnetization of the domain wall. Numer¬ 
ical simulations of the coupled equations of motion for the 
staggered field and the magnetization field confirmed that 
a spatially inhomogeneous magnetic field can act as a po¬ 
tential well for the domain wall. Finally, we have discussed 
how the intrinsic magnetization of antiferromagnetic tex¬ 
tures, which for one-dimensional domain walls is not larger 
than the spin on one atomic site, can be experimentally 
exploited in 2D and 3D. In such higher-dimensional real 
systems the intrinsic magnetization accumulates in the per¬ 


pendicular directions and the total spin can, therefore, be 
of appreciable size and may be detectable. 
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Appendix A: Energy functional for D > 1 


In this Appendix, we expand our calculation of the free 
energy functional of AFMs to 2 and 3 dimensions to disclose 
the form of the parity-breaking term in higher dimensions. 
For D = 2, we use the centered rectangular unit cell, with 
two sublattices within each unit cell. Starting with Eq. (1), 
we now assume that ol and (3 are two-dimensional vectors 
and that the coordinate pair {i, j} unambiguously defines 
all antiferromagnetic unit cells. Next, we define 


+ S 


/3 


“ 25 

Sl’^' - 


25 

+ li,,.). 

- kj) 


(Ala) 

(Alb) 

(Ale) 

(Aid) 


where we must take into account the equivalence of inter¬ 
changing \i —1^, such as in the one-dimensional deriva¬ 
tion. The Heisenberg Hamiltonian can be written as a sum 
over antiferromagnetic unit cells in the perpendicular i and 
j directions 


Ar-i,Ar-i 

H2d = JS^ - li,i)[(mi,j + lij) 

hj 

+ 1^+1,i) + 

N,N 

hj 


where we have disregarded a small energy contribution 
from the edge spins like in Sec. HA. Eq. (A2) is a sum over 
the Nn = 4 nearest-neighbor exchange couplings and the 
anisotropy energies for each antiferromagnetic unit cell. We 
use the identities — 

and {lijuii+ij - iriijli+ij) = - uiij) - 
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— \i^j) etc. to rewrite Eq. (A2) to 

AT, AT 

if2D=iV„J52^K.-l2,) 

2 N- 1 ,N -1 

ij 

+(li+i,i+i - - (mj+ij - mij)2 

-(mjj+i - mij)2 - (uii+ij+i - mijf] 

N-l,N-l 

+ JS^ E 

id 

“1“ + l “1“ 

AT, AT 

+ (A3) 

Ai 


To make the transition to the continuum limit, we define 
the derivatives in the linear approximation 

lim - hj) I [ J(l) Ai]dV , (A4a) 

hj 

^^lim ^ - hj) I [ J(l) A,-]dV , (A4b) 

^■>3 

lim ^ [[Ji\)A, + J{l)Aj]dV, 

^■>3 

(A4c) 

where J7(l) is the Jacobian matrix of the vector field 1, 
^i(j) is a vector between unit cells in the i{j) direction, and 
V is the volume of the unit cell. For the centered squared 
unit cell, \Ai\ = |Aj| = A and E = A^. We define similar 
derivatives as in Eqs. (A4) for the magnetization field m. 


The procedure is analogous when including a third di¬ 
mension, e.g., for a body-centered cubic unit cell, repeat¬ 
ing the above calculation with {i,j} ^ Apart 

from a constant contribution, the resulting free energy 
density for ATMs in dimensions D > defined here as 
H2{3)d = f(dV/V)'H2(3)D, is given by 

^2(3)D = JS'^Nn{2ui^ + - (5im)2] 

+ ^ AiAj{din ■ djU — dim ■ djm) 

+ i y] Aj(m dil-l- 9jm)} 

i 

(A5) 

where we may define i and j to run over perpendicular 
directions {x^y^z}. The sum over first order derivatives 
arises from the relation J^(1)A = A^9j(l), where A = 

{ Af, Aj, A/j}. 

By considering squared or cubic lattices, A = 2d/\fD 
and d is the nearest-neighbor distance. We can express the 
free energy density in the exchange approximation, |i^| ^ 
|J|, as 


a 2 I ^ 


^ 2 ( 3 )D = -m , ^ 


y](5in)2 + - '^diU-djU 


^7^3 


+Lj2{m-din) - , 


(A6) 


where a = 4Ar„JS'2, A = NnA^JS^/2, L = N^AJS^, 
= 2KS^, and Nn is the number of nearest neighbors. 
In antiferromagnetic materials in which the exchange en¬ 
ergy is anisotropic due to, e.g., more complicated unit cells, 
Eq. (A6) can still be used, although in this case a. A, and 
L must be treated as tensors. 
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